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Due to the good detection limits of graphite fur- 
nace atomic absorption spectrometry (GFAAS), 
situations frequently arise when an analyte is deter- 
mined in the presence of a large excess of concomi- 
tants. It is in these cases when the improvements 
made in furnace design, heating rate optimization, 
Zeeman background correction, platform atomiza- 
tion and application of matrix modifiers are still not 
sufficient to provide the analytical chemist with re- 
liable data. 

On the other hand, large gains in mechanistic 
insight have not yet yielded a model that can suc- 
cessfully cope with atomization characteristics in 
samples with a large matrix/analyte ratio. Typical 
accounts of the current state are found in [1] for 
explicit differential equations formulations and in 
[2] for the Monte Carlo approach, but there is no 
real difference in the predictive abilities of the two. 
Nevertheless, it is clear that a wide variation of 
peak shapes are observed and qualitatively moni- 
tored by everyone practicing GFAAS. 

In this work we describe a characterization of 
peak shapes by various parameters, such as appear- 
ance time, mean, mode, peak width, area/height, 
skewness and kurtosis. This is done by attempting 
to alter the peak shape systematically through de- 
liberate addition of concomitant elements in a frac- 
tional factorial experiment [3,4]. The observed 
shape parameters are then linked to the characteris- 
tic mass after discarding those that vary with ana- 
lyte mass. 



Experimental 

A Perkin-Elmer 50Q0 with HGA 500 and AS 40 
autosampler was used throughout. The computa- 
tions were run on a Perkin-Elmer 7500 data station 
and on an IBM PC-XT. All data given are Zeeman 
background corrected and all samples are atomized 
from a platform. Two different platforms were 
used. One was the vendor-supplied pyrolytic 
graphite platform, and the other one was produced 
from the same type by physical vapor deposition of 
TaC [5], A general account of the analytical prop- 
erties of this TaC-coated platform for ETAAS can 
be found elsewhere [6]. All reagents were of ana- 
lytical grade. Doubly distilled water and subboiled 
nitric acid were employed. 

Sn, Se, and Rh served as model analytes to study 
the behavior of elements different in both chem- 
istry and volatility. Ten micrograms of each palla- 
dium and magnesium as nitrates were used as 
matrix modifiers for Se, and 200 u-g ammonium di- 
hydrogen phosphate was used for Sn. First, tem- 
perature programs were optimized for standard 
solutions and not altered for any of the subse- 
quently described experiments. Generally, atomiza- 
tion temperatures from the TaC-coated platform 
tended to be higher than for the ordinary total py- 
rolytic graphite platform [5]. 

Calibration curves were run to establish the lin- 
ear range, but — more importantly — the depen- 
dence of peak shape parameters on analyte mass. 
From this the concentration levels were chosen for 
the interference experiments to be well within the 
linear range. Cr, Fe, Na, Zn, Mg and Al were cho- 
sen as potential interferents and added as nitrates in 
a 29000 (Sn) to 40000 (Se) fold molar excess over 
the analyte. A fractional factorial design in these 
interferents was set up which was replicated once, 
giving a total of 64 profiles [4]. These were 
recorded in random order with 2 atomizations of 
the standard solution and 2 atomizations of a blank 
solution before and after the run of 64. 



Results and Discussion 

Peak shape parameters can serve as a reduced 
description of the entire atomization curve and 
should be defined to carry information on the 
(largely) unknown fundamental mechanisms in 
GFAAS. In this work it was assumed that— if not 
the mechanisms— at least kinetic and thermody- 
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namic parameters governing the processes would 
be altered by the presence of interferents and result 
in changes of peak shape. Consequently, the ob- 
served atomization curve could be regarded as a 
faithful representation of the free atom population 
in the detection volume and should be accessible to 
a nonparametric description. The peak shape 
parameters used in this study are given in table 1, 
As these parameters were chosen to model sensitiv- 
ity in unknown samples, it is required that they be 
fairly insensitive to analyte mass. This could be 
shown to be the case for all, but the parameter 
"length" [7] which was therefore not used in fur- 



data was randomly removed for the modelling step 
and used to assess the predictive ability afterwards. 
Table 2 gives an overview of the results. The opti- 
mal number of dimensions for the PLS model 
varied between 2 and 6 with 62 profiles fitted; six 
atomization curves were used for predictions. The 
relative importance of the different peak shape 
parameters is different for each data set. Detailed 
results on this subject will be published elsewhere 
[9]. 

Table 2. Results of modelling characteristic mass from peak 
shape parameters 
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time interval from appearance time 
to drop of the signal below signif- 
icance 

width of the peak measured by 
assuming the signal to be the prob- 
ability density curve of the free 
atom population in the analysis 
volume 

time of peak maximum 

arithmetic mean time of the peak 

measure of asymmetry of peak [3] 

measure of slimness of peak [3] 

absorbance corresponding to the 
average spread of baseline noise 



Figure 1 shows the change in kurtosis with the 
type of interferent and with the resulting character- 
istic mass. Of course, in a real analytical setting 
situations may arise where the interferent is not 
identifiable. It is therefore necessary that the char- 
acteristic mass is modelled independent of the na- 
ture and mass of an interferent. This can be 
attempted using only the peak shape parameters. 

The technique employed for this modelling is 
known in the literature as partial least squares 
(PLS) regression [8]. The data from the interfer- 
ence study, supplemented by the four data on the 
standards were subjected to PLS. A subset of the 



The success of these initial experiments is obvi- 
ous. In practice the applicability of the proposed 
scheme will depend (i) on whether the majority of 
changes in peak shape can be observed in the initial 
experimental design, and (ii) on the validity of the 
implicit assumption that most interferences mani- 
fest themselves by changes of peak shape. Both 
questions cannot be answered without full insight 
into the mechanisms of atomization. In the mean- 
while, the approach can serve as a means to warn 
the analyst about the unexpected occurrence of in- 
terferences. The model for characteristic mass can 
serve to improve accuracy. This by itself should 
help to make data from ETAAS more reliable. 
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Figure 1. Variation of kurtosis with nature of analyte and 
change of characteristic mass. 
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Introduction 

Drug abuse is a serious problem in the world and 
these drugs are synthesized in clandestine laborato- 
ries. The determination of the manufacturing 
source provides important information regarding 
drug traffic [1,2]. One of the stimulants imported 
illegally, methamphetamine, is a serious social 
problem in Japan. 

Methamphetamine is synthesized from ephedrine 
or methyl benzyl ketone as the raw material. In 
Asian areas, ephedrine is used as a starting material 
and methamphetamine is prepared by catalytic re- 
duction of chloroephedrine by palladium-barium 
sulfate [3] or by reduction of ephedrine with hy- 
drogen iodide and red phosphorus [4] as shown in 
figure 1. 

The identification of catalyst and/or reagents 
based on final product analysis have been studied 
by radiochemical neutron activation analysis 
(RNAA) [2]. Palladium, barium, and iodine were 
determined as trace elements in methamphetamine 
hydrochloride. 

In this paper, these trace elements were deter- 
mined directly by inductively coupled plasma-mass 
spectrometry (ICP-MS) and the results are com- 
pared with those of neutron activation analysis 
(NAA). 



Experimental 

Reagents 

Methamphetamine hydrochloride was prepared 
by Emde's method [3] or Nagai's method [4] as 
shown in figure 1. 

Procedure 

For the ICP-MS, samples were dissolved in wa- 
ter in a concentration of 0.1-1%, and then the ana- 
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